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INTRODUCTION 


_ The data reduction utility METDAT enables interpretation of raw 
meteorological data from the USDA Forest Service EMCOT Weather 
Station portable towers (Ekblad et al. 1990) into statistical averages for 
input into FSCBG or graphical presentation packages. The structure of the 
data logger data files are incorporated into the program to read and 
interpret the raw data within a simple user interface. Meteorological 
Statistics are then generated for the time intervals desired. The 
accompanying 3-1/2 inch diskette contains the program for operation on a 
personal computer. 
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METDAT OPERATION 


METDAT may be initiated with the command 
METDAT <cr> 


followed by a carriage return < cr > or enter keystroke. The program 
must be in the current subdirectory, or accessed through the DOS PATH 
statement. Alternately, METDAT may be initiated with a filename on the 
command line 


METDAT metrawfile < cr > 
where metrawfile is the name of a comma-delineated ASCII file containing 
the data logger recorded data (as from Thompson 1991). By convention 
this file extension is always DAT. 
METDAT is configured for operation on CRT color monitors; thus, 
LCD gray-scale monitors may fail to represent some menus adequately. 
To overcome this effect the user may initiate METDAT in a black-and- 


white mode by including a separate character string ("/g") on the command 
line to indicate operation on a gray-scale monitor: 


METDAT /g <cr> 
METDAT metrawfile /g <cr> 
All operations will then commence in black and white only. 
In succession, METDAT makes the following screen entries: 
Enter name of met raw data file 
will be asked if metrawfile is not included on the command line. This file 
should be in the same subdirectory as METDAT. If an incorrect file name 
is entered, the error message 


Unable to open metrawfile. DAT 


will appear. Subdirectory contents should be examined to select a correct 
file name. 
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Checking metrawfile. DAT 


will be shown as the program reads through the data file (recorded times in 
HoursMinutes_ will be displayed to indicate that METDAT is reading the 
data file). 


Time limits are between xxxx and yyyy 


will indicate what the actual time limits are in the data file) METDAT will 
generate ten-minute averages of meteorological data, beginning with a 
complete time recorded by the data logger (a time that includes 
temperature, relative humidity and wind data). 


Some of the met raw data files from Thompson (1991) appear to be 
corrupted at the beginning of the file. METDAT attempts to overlook this 
data by checking for consistency from the end of the file backwards to the 
front. 


Enter time to begin averaging 


will be asked to select the first ten-minute interval for averaging. The time 
selected is the time at the beginning of the interval, and must be within the 
time limits given previously. 


Enter number of averages to compute 


will be asked to indicate how many ten-minute averages should be 
performed. The averages are subsequently displayed on the screen, and 
read into a comma-delineated ASCII file (with the extension AVE) for 
transfer to statistical and graphical packages such as Lotus 1-2-3 or 
Harvard Graphics. The screen will display the first ten averages. Once the 
message 


zZzzz Entries written to metrawfile.A VE 


appears, the user may examine this file, on the screen, with the aid of the 
cursor up or down arrow keys. A carriage return < cr > or enter will exit 
the program. 


The screen presents the column results (duplicated in the created 
ASCII file metrawfile.AVE but without the units) as follows: 
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Time The HoursMinutes designation 
from the data logger for the time 
at the beginning of the ten-minute 
averaging interval 


Wind Speed (in meters/second) The average horizontal wind 
speed in each ten-minute interval 
(20.0 feet above the ground) 


Wind Direction (in degrees) The average wind direction in 
each ten-minute interval (20.0 feet 
above the ground) 


Sigma A (in radians) The azimuth standard deviation of 
the horizontal wind speed 


Sigma E (in radians) The elevation standard deviation 
of the vertical wind speed 


Temperature Bottom (in deg C) The average temperature at the 
lower temperature sensor (4.0 feet 
above the ground) 


Temperature Top (in deg C) The average temperature at the 
upper temperature sensor (20.0 
feet above the ground) 


Relative Humidity (in percent) The average relative humidity 
(20.0 feet above the ground) 


In all cases Sigma E has been set equal to one-third of Sigma A. 
This approach is a compromise for the inherent unreliability of the 
propeller anemometer to detect vertical wind speed (Bowers 1991). 
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STATISTICS 


Mean horizontal wind speed U is determined from the equation 
U=1>DU; 
N 


where the index i denotes each data logger entry for wind speed in the ten- 
minute interval, and N denotes the total number of data points considered 
(typically N = 600 for data recorded every second). 


Mean vertical wind speed W is determined by the equation 
w=4dw; 
N 
where the data logger entries for W-component are used. 


Mean wind direction © is found from the unit vector approach 
(Haugen 1963) by defining 


X =1Y cos @; 
N 
Y =1Y sin @; 
N 
to obtain 
© = arctan | x] 
xX 


where, internally, METDAT uses the ATAN2 intrinsic function defined on 
the interval -180 degrees to 180 degrees . This procedure eliminates any 
erroneous behavior in the wind direction because of cycling of the 
direction angle every 360 degrees. It is anticipated that the meteorological 
tower will be set in the field so that North corresponds to © = 0 degrees. 
Thus, © determined here is the wind direction from North, with a 
positive value indicating clockwise. Other tower positions will change the 
zero position of the wind direction accordingly. 





















2 ITZITATS 


te 5 ole aT toes e. 
noitsups st mot bentansieb af U besqe barw lemoxhet mes 
. ww aie 


13) cli ai bssqe baw 103 Yuis 19920! } mn fost extonsh { robot sett & 
hsp mo) 2Aniet feb ° iO wd vite Lalo) nt B 1 BST: 123 vi bras avy 4 ! 


(bacsse vers betroos ate 10? ee 
noitenpe ont yd barierrsiab 7 be ge horw ae 


iW Sh = fi 


' - - - " 
beau sis tratoumoo-W it esitine t9ggol ateb of. ore 
pee: 
me Toloev tin ori ott bavot 2: © aolssedb baiw ns ; 


-— 


; Lex) 
i9 209 or 4 


Eleee, 


no benitsb notional sianinali SMATA ont sago TAQTSM , lees mn un 9 
(ne asioaimils swbevor 2idT . esoxgsb 081 of preatr 
ett? to prllsyo to eausosd. poiszerib bniw art of rived ad ‘sian . 
Isvigolmostain od jer hateqioitas af I pennant nav se 
2osigth O = @ ot ehaogeonco dnoM tari poe 2 od 

6 diw .dnold men? aoitoodb boiw od} ai on: S198 ‘etanaae a 
att gmat iw ancisnog 12 nfo gait 


Om det 9 
5 


“ 

i 4 Sit): 
Pah a ’ 
easel a on ag 
aa a TV : 


sud DAT 





28 “7 
i. .. ra 7 


LP EP, g Oh 


Azimuth standard deviation 6, is then found from the equation 
(Haugen 1963) as 


ee 1/2 
SA = Aas (1-x2-¥| | 
while elevation standard deviation Og is found from the equation 
op=to 
Be 3 A 
to avoid propeller anemometer signal difficulties at low wind speeds. 


Finally, average temperature T and relative humidity R are found 
from the equations 


T=1>T; 
Ne : 
and 
R=1ER; 
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